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Abstract 

A modeling of the integral-type Stirling refrigerator (ITSR) using a system dynamics approach is developed in the 
present study. The system dynamics models are derived to describe the input/output relation of each component of an 
ITSR. Connecting the equivalent circuits of the components together, we obtain a flow network diagram and derive a 
transfer function to represent system dynamics behavior of an ITSR. The performance analysis of an ITSR is then 
easily carried out by use of the transfer function with sinusoidal signal assumption. The performance calculation for a 
test refrigerator has shown that the analysis is accurate. An empirical correlation for a correction coefficient r in the 
expansion space was also derived for better results. The present study has demonstrated how a linear system theory or 
system dynamics approach can be applied to the design of an ITSR. © 2000 Elsevier Science Ltd and HR. All rights 
reserved. 
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Moderation d’un refrigerateur Stirling de type integre a 
l’aide d’une approche par dynamique du systeme 


Resume 

Dans cette etude, on a effectue la modelisation d’un refrigerateur Stirling de type integre (ITSR) a I’aide d’une 
approche par dynamique du systeme. Les modeles de la dynamique du systeme sont developpes afin de decrire la relation 
entree/sortie de chaque composant d’un ITSR. En raccordant les circuits de composants, on obtient un diagramme pour 
I’ecoulement en reseau avant de deriver line fonction de transfert qui represente le comportement dynamique du systeme 
d’un ITSR. L'analyse de la performance d’un ITSR devient alors facile grace a I’utilisation de la fonction de transfert 
basee sur un signal sinusoidal. Le calcul de la performance d’un refrigerateur d’essai a ensuite montre que cette analyse est 
precise. Une correlation empirique pour un coefficient de correction r dans la chambre d’expansion a ete derivee afin 
d’obtenir cle meilleurs resultats. Cette etude a montre qu’une theorie d’un systeme lineaire ou une approche par dynamique 
du systeme peut etre appliquee ci la conception d’un ITSR. © 2000 Elsevier Science Ltd and HR. All rights reserved. 

Mots cles: Cryorefrigerateur ; Stirling ; Conception ; Modelisation 
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Nomenclature 

Greek letters 



£ 

regenerator porosity, dimensionless 

^dp 

gain of G dp (s) 

T 

time constant (s) 

^ep 

gain of G ep (s) 

CO 

angular frequency ( = 2 nf) (rad s" 1 ) 

A 

area (m 2 ) 

P 

density (kg m -3 ) 

C 

flow capacitance (m 3 kg kJ -1 ) 

0 

displacer leading angle to piston, deg 


specific heat at constant pressure (kJ kg~' 




K- 1 ) 

Subscripts 

C v 

specific heat at constant volume (kJ kg^ 1 

c 

compression chamber; characteristic impe¬ 


K- 1 ) 


dance; cold end 

h 

convection heat transfer coefficient (W nr 2 

ce 

cold end 


K- 1 ) 

d 

displacer 

d 

diameter (m) 

e 

expansion space 

f 

frequency (Hz) 

f 

gas; fluid; flow 

k 

thermal conductivity (W m -1 K~*) 

g 

gas 

L 

length (m) 

h 

hot end; hydraulic 

m 

mass flow rate (kg s“ x ) 

H 

hot end of the regenerator 

R 

gas constant (kJ kg^ 1 K _1 ) 

i 

inlet 

Re 

Reynolds number, dimensionless 

L 

cold end 

P 

pressure, Pa (N m~ 2 ) 

m 

mean; regenerator matrix 

T c h 

charge pressure, Pa (N m~ 2 ) 

o 

outlet; equilibrium point; amplitude 

Q net 

net cooling capacity (W) 

P 

piston 

S 

Laplace variable, complex 

r 

regenerator 

T 

temperature (K) 

s 

solid material 

t 

time (s) 

t 

gas passage 

V 

volume (m 3 ) 

w 

warm space; cylinder wall 

X 

distance (m) 



X 

displacement (m) 

Superscripts 




perturbed value 


1. Introduction 

An integral-type Stirling refrigerator (ITSR) com¬ 
prises a compressor and a cold finger that are connected 
by a gas passage. The cold finger contains a moving 
regenerator/displacer, a cylinder, and an expansion 
space as shown in Fig. 1. Both the piston and the dis¬ 
placer motions are driven by a crank-shaft mechanism 
with a fixed phase angle between them. An electric 
motor drives the piston and hence the displacer such 
that the refrigerator operates approximately in a Stirling 
cycle. The physical processes involved in an ITSR 
includes thermodynamics, heat transfer of porous med¬ 
ium in the regenerator, gas transport in connecting 
tubes, tribology for the friction of piston and displacer 
motion. In addition, all the processes are at unsteady or 
cyclically-steady state. The thermal analysis of an ITSR 
is thus quite complicated. 

Considerable work has been done by many research¬ 
ers on the cycle analysis of the Stirling machine. The 
analysis basically utilizes the principle of thermodynamics 
and the conservation laws for mass, momentum and 
energy. The analysis is classified into two categories by 


Organ [1], engineering thermodynamics approach and 
flow-field approach. In either approach, a numerical 
calculation using the finite difference method is sometimes 
employed to calculate the cycle performance and the form 
of the solution is mostly for a single operating point in the 
performance envelope of a specific machine. A perfor¬ 
mance map for a specific machine is hardly obtained for 
the selection of an optimum operating condition. 

Urieli and Berchowitz [2] treated a free-piston Stirling 
machine as a generalized single degree of freedom sys¬ 
tem (damped mass-spring system in mechanics) and 
derived the general relations for piston-displacer ampli¬ 
tude ratio and phase angle in terms of operating fre¬ 
quency. In order to carry out the cycle performance 
calculation, Urieli and Berchowitz [2] further used the 
closed-form solution of the isothermal analysis for the 
working gas pressure as a function of compression and 
expansion volume. The relations of piston-displacer 
amplitude ratio and phase angle for a specific machine 
were derived and used in the cycle performance calcula¬ 
tion. Since the isothermal model is an ideal model, the 
analysis of Urieli and Berchowitz [2] for a free-piston 
Stirling machine will be inaccurate. 
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compression 

chamber 



space 


COMPRESSOR COLD FINGER 

Fig. 1. Schematic diagram of an integral-type Stirling refrig¬ 
erator. 

Fig. 1. Schema d’un refrigerateur Stirling de type integre. 


ITSR is directly driven by a motor through a mechan¬ 
ism such as a crank-shaft mechanism. The displacer has 
a fixed stroke and a fixed phase angle with respect to the 
piston motion. For a STFSR, the displacer is driven by 
the fluid force or gaseous pressure which is indirectly 
induced by the piston motion. The stroke is thus vari¬ 
able and depends upon the operating conditions. 

In modeling the ITSR by using a system dynamics 
approach, the piston displacement is the system input. 
The integral of the gas pressure in the expansion space 
with respect to the displacer motion is the gross refrig¬ 
erating power that the refrigerator can obtain. The 
expansion-space pressure is, thus, the system output. 
That is, the refrigerator can be treated as a dynamical 
system with the piston displacement X p as the system 
input and the gas pressure in the expansion space p e as 
the system output. A transfer function representing the 
performance of the refrigerator can then be derived, 
from which the system performance can be easily eval¬ 
uated at various operating points. 


For the analysis of an ITSR, a nodal model based on 
numerical computation scheme has been developed [3]. 
Several commercial software packages are available. 
However, the computation suffers easily from numerical 
problem and the accuracy in design is not known. In the 
present study, we intend to derive a model for the per¬ 
formance analysis of an ITSR that is different from the 
conventional approach in concept. A closed-form solu¬ 
tion is obtained in transfer-function form and from 
which the performance of an ITSR is easily calculated. 

Since an electric motor drives the piston and the dis¬ 
placer in sinusoidal motion approximately, a cyclically- 
steady state operation very close to sinusoidal is 
obtained in an ITSR. Every component of the refrig¬ 
erator is thus subject to a sinusoidal signal and gen¬ 
erates a sinusoidal output just like the frequency 
response of a dynamical system, from the point of view 
of system dynamics. A dynamics model can then be 
derived to describe the input/output relation of each 
component using the governing equations based on the 
conservation laws for mass, momentum and energy and 
the principle of thermodynamics. A linear system 
dynamics model can be further derived in conjunction 
with linearization and approximations. A system dynam¬ 
ics model and an equivalent system block diagram or 
equivalent circuit can then be used to represent the 
dynamics behavior of each component. Connecting the 
analogous circuits of all the components together will lead 
to a network representing the operation of the refrigerator. 

Huang and Lu [4] derived two transfer functions for 
the split-type free-displacer Stirling refrigerator (STFSR). 
In the present study, we derived a transfer-function 
model of the ITSR for the system performance analysis. 

The major difference between ITSR and STFSR is the 
driving mechanism of the displacer. The displacer of an 


2. Dynamics model of components 

Assuming that the working fluid is an ideal gas, we 
obtain the dynamics models of each component of an 
ITSR which are summarized as follows. 


2.1. Compression chamber 


Assuming an isothermal compression with a constant 
gas temperature T c and no gas bypass between the pis¬ 
ton and the cylinder, we obtain a dynamics model: 


m p (s) 


sp 0 A 

~RT C 


-IpW = m c (s) 4 


Pc(s) 

1 /(sC c ) 


( 1 ) 


where C c = V co /RT C is the capacitance of the compres¬ 
sion chamber. Here, m v (s) acts as a current source 
representing the available or gross mass flowrate gener¬ 
ated by the piston motion; m c is the mass flowrate out of 
the compressor; p c is the gas pressure inside the com¬ 
pression chamber; V c (t) = V co — A p X p (t); V co is the 
volume of the compression chamber for the piston at 
the equilibrium position; X p (t) is the piston displace¬ 
ment measured from the midpoint of the piston motion 
toward the top dead end; A p is the cross-section area of 
the piston. 

2.2. Gas passage 

The gas passage links the compressor and the cold 
finger. Following the same treatment of Huang and Lu 
[4] for a connecting tube (assuming an isothermal pro¬ 
cess with a constant temperature T t ), we obtain a 1-D 
dynamics model of the passage: 
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riho(s) 


cosh(r t L t ) — sinh(r t L t ) 

_ Pto(s) _ 


—Z ct sinh(r,L t ) cosh(r t L t ) 


»2tiO) 

Pds) 


( 2 ) 


where T t = \/Z t Y t and Z ct = s/Zf]Y[ are the propaga¬ 
tion constant and the characteristic impedance of the 
passage, respectively; Z t = R Ft + .sL Ft = series impe¬ 
dance and Y t = sC Ft = shunt admittance; 

r t = yjsC\ fR\ v + sLp t ); Z ct = ,, 1 (3) 


where C F t, Lpt and Apt are the flow capacitance, flow 
inductance, and flow resistance per unit passage length, 
respectively, which are defined as 


tilr oO) 
Pvo(s) 


— R(s) x 


m ri (» 
Put*) 


where 


m = 


Rmm (s) 
Rpi]] (p) 


Rmp (•S’) 

RppO 5 ) 


(7) 


( 8 ) 


— 7?pp(^) — COSh[r r (5)L r ] 
sinh[r r (s)L r ] 

" p<l) - Z tr (s) 

/'Jpmfv) = —Z cr (.v) sinh[r ,f.v) /.,] 


(9) 


where T r = \fZfT x and Z cr = sfZfJT x are the propaga¬ 
tion constant and the characteristic impedance of the 
regenerator, respectively; Z r = R Fi + sL Fi = series impe¬ 
dance and Y r = sCFr = shunt admittance; 


CFt = yyyr; — —; R Ft = — 

RT t A t A t 


(4) 


r r = yj ,S'C| Tr( R\ \ + sL F i )', Z cr = 


r r 


■VC|Ti 


( 10 ) 


where 

K = 0.1556(p t w max dt/At) -0 - 201 (vrmax/4), (5) 

A t is the cross-sectional area of the passage; w max is the 
peak velocity of the oscillating flow in the passage; dt is 
the hydraulic diameter of the passage, k is considered to 
be a constant during the modeling, but should be 
adjusted by numerical iteration during the computation 
[4] to give a correct value for the corresponding mass 
flow H’ ma x- 

2.3. Warm space 


where fj| r is the regenerator flow capacitance due to 
pressure change and time responses of gas and matrix: 


^ , i 1 4 * Ar/[ T gr(l T sr sr )] 

FTr “ Fr y + W[r gr (l+.vt sr )] 


( 11 ) 


where y = C p /C v ; C Fr and Cj r are the flow capacitance 
per unit length due to pressure and temperature changes, 
respectively; L Fr and R Fr are flow inductance, and flow 
resistance per unit length, respectively. They are defined as 


,, 4 f r PrmA'v T 1 „ USV 

Cpr - RT, m ■ Cti ~ R7%.' F ‘ “ A s , ’ Rf ’ - A„ 


( 12 ) 


The gas in the warm space is assumed to be well- 
mixed at an uniform temperature. Similar to the com¬ 
pression chamber, a linearly perturbed model can be 
derived from mass continuity Eq. (4): 


ih w fs) 


spo A (jw 
RT W 


Xd(s) + m WO C ?)4 


Pw(s) 

l/(.s'C w ) 


( 6 ) 


where C w = V wo /(RT w ). The first term in the right-hand 
side of Eq. (6) is the mass flowrate induced by the dis¬ 
placer motion, which is analogous to an electric current 
source; the third term represents the flow capacitance 
effect of the warm space. 

2.4. Regenerator 

The regenerator is an energy-storage element made 
from packed wire screens. Assuming 1-D flow, no axial 
conduction and constant properties, Huang and Lu [4,5] 
obtained the dynamics model of the regenerator: 


where 

ci = ci + (/3/(/h)£Reh; a = 175/(2e^); f) = 1.6/(2e 2 dh); 
dp — ed m /( 1 — e); h = 0.33(kf/t/h)Re{|' 67 

_ PrmC i 'of \ _ _ Ps(l £') I j)fs 

Tgr _ RT rm hA m ; Tsr ~ hA^j 


Aut is the surface area of the regenerator = 
4To(l —s)/d m ; Reh is the Reynolds number based on 
hydraulic diameter <4; /tf is the thermal conductivity of 
the gas. 

Eq. (7) represents the dynamics model of the regen¬ 
erator in terms of the pressure and mass flowrate that 
are derived from the continuity and momentum equa¬ 
tions. The gas temperature at the cold side of the 
regenerator is derived as [4] 
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TfALr, s) = - 


CFr Tg r ( 1 + ST sr ) 
Cjr T sr 


Pr oO) 


Y r gr (l+5t sr )r r (5) 
Cjr sr ST sinh[r r L r ] 


x [m ro (s) cosh(r r L r ) — /» n (.v)]. 


(13) 



2.5. Displacer 


COMPRESSION GAS WARM REGENERATOR EXPANSION 

CHAMBER PASSAGE SPACE SPACE 


The displacer of an ITSR is driven by a crank-shaft 
mechanism (Fig. 1). The time functions for the displacer 
and piston motion are approximately sinusoidal. 
Therefore, 


Td(0 = Td 0 sin(ft>? + 0) 

Tp (/) = X po si w{(oi) 


where X do and X po are the amplitudes of the displacer 
and the piston, respectively; 6 is the phase lead of the 
displacer. Thus, a transfer function G dp (s) for the dis¬ 
placer motion can be derived as 


G dp (s) = 


Xd(s) 

*p(s) 


X d Q ,n 

-e 

*po 


(14) 


Fig. 2. Linear network of integral-type Stirling refrigerator. 
Fig. 2. Reseau lineaire du refrigerateur Stirling de type integre. 

The lower part of Fig. 2 is the fluid network in which 
mass flow is analogous to current and pressure is ana¬ 
logous to voltage. The fluid network depicts the gaseous 
flow propagation from the compressor to the expansion 
space. The upper part of Fig. 2 represents the displacer 
model which responds the piston motion by a phase 
angle 6 and transfers the displacer displacement signal 
X d to the fluid network. Since the displacer motion is 
fixed and follows the piston motion by a phase angle, we 
can solve the fluid network to find the gas flowrate and 
pressure at each component. 


2.6. Expansion space 

The expansion space is usually an empty space for the 
gas to expand and absorb heat from the environment. 
Huang and Lu [4] derived a dynamics model from the 
mass continuity equation. To take into account the gas 
bypass from the displacer clearance, we modify the flow 
capacitance C e (= V eo /RT ee ) with a correction coeffi¬ 
cient r and obtain the dynamics equation: 


4. Transfer function of integral-type Stirling refrigerator 

We draw only block diagram in Fig. 2 for the gas 
passage and the regenerator since their dynamics models 
belong to distributed-parameter type. The transfer 
function of the ITSR is derived by combining the above 
equations: 

Ge P (s) = = G 9 (s ) + 6- |0 (.v)G dp (.v) (16) 


m c (s) 


Pe(s) 

1 /(rsC e ) 


SPo 1 clc 

ET ce 


X d (s) 


(15) 


The second term in the right-hand side of Eq. (15) is 
the mass flowrate induced by the displacer motion, 
which is analogous to a current source. The correction 
coefficient r represents the gas bypass from the displacer 
clearance. The amount of gas bypass is related to the 
seal design, manufacturing process, material used and 
operating condition. The relation of r can be determined 
empirically, r = 1 represents no gas bypass in the 
expansion space. 


where 

C M 4f c i(.?) X do e 

G ^ s) -x^)-x;j 


GioC?) — 


Gg(s) = 


G 6 (s) - G 7 (s)E md (s) 
G 8 (s) + G 7 (s)E mc (s) 


1 


G 8 (s) + G 7 (s)E me (s) 

G 8 (s) = -Gfs)R mp (s) + G 5 (s)R mm (s ) 
G 7 (s) = Gfs)R pp (s) - G s (s)R pm (s ) 


(14) 

(17) 

(18) 

(19) 

( 20 ) 


3. Analogous network of integral-type Stirling refrigerator 

Combining the equivalent circuits or block diagrams 
of all the components according to the process of ITSR, 
we obtain the analogous network as shown in Fig. 2. 


G(fs) = CE(s)W nui (s) 

(21) 

G 5 (s) = G 4 (s) - Gfs) W mp (s) 

(22) 

G 4 (s) = -G](s)f mp (s) + G 2 (.s')/ mm (.s') 

(23) 
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(his) = 

CiCs)/pp(.s) - G 2 (s)I pm (s) 

(24) 

G 2 (s) = 

Cinc(s) _ 
Cmp(^) 

Giis) = 

1 

CmpC?) 

(25) 

C m p(s) : 

sp 0 Ap 
" KT C ’ 

C m c (A 

ii 

C-0 

(26) 

fFmp is) 

s Fwo 

RT W 

; w md (s) 

sp 0 A d w 

RT W 

(27) 

Eme(^) = 

_ rsV e o _ 
RT c e ’ 

'mdCs) = 

sp 0 A de 

RTce 

(28) 

7mm Of) = 

- 7p P (.s') — 

cosh[r t (.s)L t ]; 

(29) 


Peo — | Cep (jt'P 11 A p ( j( 0 ) | — a ep A p0 
do = \G dv (jco)\\X v (ja>)\ — a dp X po 


(32) 


The maximum cooling capacity (7 m: , y thus can be com¬ 
puted simply from the gains of the frequency response 
functions G ep (jco) and G dp (jco), i.e. a ep and a dp , and the 
phase angles of X d (t) and p e (t ), i.e. fa and fa: 


gmax =/ 


2n 


PedV e =firp eo A de X do sm(fa - fa) 

= fita cp a dp X 2 po A Lk sin(0 d - fa) 


(33) 


where A de is the cross-section area of the displacer. 


/mp(.s') - - smh ^ Ll] ; / pm (. y ) = -Z ct sinh[r t (s)L t ] 

(30) 

It is seen that the system input of an ITSR is the piston 
displacement X p of the compressor. The system output 
is the expansion space pressure p e . The integral-type Stir¬ 
ling refrigerator thus belongs to a single-input-single- 
output (SISO) dynamical system. 


5.2. Net cooling capacity 

The net cooling capacity of an ITSR can be evaluated 
by subtracting the heat losses from the maximum avail¬ 
able cooling capacity. There are three types of heat losses: 
namely, regenerator heat conduction loss g con d, regen¬ 
erator enthalpy flow loss g en th> and displacer shuttle 
heat loss g s hutt- Thus the net cooling capacity is 

Qh = gmax gshutt gcond genth (34) 


5. Performance calculation 

5.1. Maximum available cooling capacity 

The maximum available cooling capacity of an ITSR 
can be evaluated by integrating the expansion space 
pressure p e (t)[= Pe + p e {t )] with respect to volume V e (t), 
where V e (t) is related to the displacer displacement 
X d (t). For a given piston and displacer motion, X p (t) 
and X d (t), according to the crank-shaft mechanism 
design, p e (t) can be obtained from the inverse of the 
transfer function G ep (s), Eq. (16). The computation 
procedure is however quite complicated. A simplifica¬ 
tion is needed. Since the piston is driven by a crank¬ 
shaft mechanism, the piston and the displacer motion as 
well as the induced pressure and mass flow waves inside 
the refrigerator approaches sinusoidal with the time 
functions of X d (t) and p e (t) expressed as 


where 


gcond = /c e ff(7T.D;/4) 


T h -T l 

L r 


2enth —/(D /^ro( 0 Cp Tf r (L r 5 t)dt 


— fnC p m ro | 


Tfr(L r ) 


COS(0 m - fa) 


(35) 


(36) 


where is the effective thermal conductivity of the 
regenerator; D r is the regenerator diameter; L r is the 
regenerator length; 7 h and T f are the temperature at 
the hot and cold ends of the regenerator; <p m and fa are, 
respectively, the phases of m IO (t) and Tf r (L r , t). 
T[ r (L r , t) follows Eq. (13). The displacer/regenerator 
shuttle heat loss is evaluated from the relation of Baik 
and Chang [6]: 

gshutt = Tl[) ^ Sc 2:rf k d p d C pd T l0 sirup (37) 


Pe(t) = Peo sin(6» + fa) 


(31) 


x d (t) = X do sin(6» + fa) 


where p eo and X do are the amplitudes of the expansion 
space pressure and the displacer, respectively; fa and fa 
are the phase lead angles of p e (t) and X d (t) with respect 
to the piston motion. The amplitude of the expansion 
space pressure can then be expressed as 


where 


Tio = A T d S e 


a 2 


a\ + a 2 \ 2 , fa\ + a 2 , x 2 

+ ( -7=-h Cl\Ch 


V2 


V2 


4> — tan 


-l 


(a\ + a 2 \ 

-7=-hUl«2 

V2 _ 

a\ + a 2 

V ~7T J 


(38) 


(39) 



638 


B.J. Hucmg, H.Y. Chen / International Journal of Refrigeration 23 (2000) 632-641 






(40) 


where D e is the displacer outside diameter; S e is the 
displacer stroke; Aq is the displacer thermal conductivity; 
k w is the cylinder wall thermal conductivity; pd is the 
displacer density; C p d is the displacer heat capacity; A I’d 
is the axial temperature gradient of the displacer; ad and 
a w are the thermal diffusivity of the displacer and the 
cylinder wall; h is the heat transfer coefficient. 

The gas bypass loss at the clearance between the dis¬ 
placer and the cylinder wall is related to the seal design, 
manufacturing process, material used, and operating 
condition. A correction coefficient r is introduced in the 
expansion space model, Eq. (15), to account for this 
elfect. 


5.3. Refrigerator COP 


The COP of the integral-type Stirling refrigerator can 
be evaluated by the relation: 

C0P = w L ~ (41) 

W comp 

where 


W, 


comp 


= fi>p c cN c 


(42) 


For simplification, the compression work computed 
from the above equation does not take into account the 
motor efficiency. 


6. System design analysis and experimental verification 

The system design analysis of an ITSR can be carried 
out using the present model. For simplification, the fol¬ 
lowing assumptions are made in the analysis: 

1. Since the mixing is severe in the compression 
space, the gas temperature in the compressor can 
be assumed uniform. The gas temperature inside 
the gas passage is also assumed uniform since it is 
very short in ITSR. Therefore, the regenerator 
inlet temperature T r \ approximates the gas tem¬ 
perature in the compression chamber T c and in 
the passage T u i.e. T ri « T x « T c . 

2. The gas temperature of the warm space is con¬ 
stant and equal to the outside wall surface tem¬ 
perature Tu and the regenerator inlet temperature 
T n . That is, T w « T n « Tu- 

3. The gas temperature in the cold space (expansion 
space) Tee is approximately equal to the outside 
wall temperature of the cold end, T\_. That is, 
T ce ~ 7l. 


4. The mean pressure p Q is equal to the charge pres¬ 
sure P ch . 

5. The wall temperatures at the two ends of the 
regenerator, T n and T m , are approximately con¬ 
stant. Thus, T w « T n ~ T t « T c : T ro ss T cc ~ T, . 


Table 1 

Design specifications of a test refrigerator 
Tableau 1 

Specifications pour la conception d’un refrigerateur d’essai 


1 . 

Compressor 

Piston diameter 

14 mm 


Piston stroke 

2 mm 


Connecting rod length 

17 mm 

2. 

Gas passage 

Length 

25 mm 


Inside diameter 

1.5 mm 

3. 

Warm space 

Volume 

0.5 c.c. 

4. 

Regenerator 

Wire screen material 

SS316 


Wire screen mesh no. 

500 


Wire screen length 

37 mm 


Wire screen diameter 

4 mm 


Number of screen disks 

740 

5. 

Displacer 

Stroke 

2.0 mm 


Connecting rod length 

20.6 mm 


Diameter 

6.15 mm 

6. 

Expansion space 

Diameter 

6.15 mm 


Dead space length 

1.5 mm 

7. 

Design-point operating conditions 
Charge pressure 

13.6 kg/cm 1 2 3 gage 


Frequency 

35 Hz 


Cold-end temperature 

80 K 


Net cooling capacity 

0.2 W 


Table 2 

Measured and calculated cooling capacity at P ch = 13.6 kg/cm 2 
gage; f = 35 Hz 


Tableau 2 

Puissance frigorifiaue mesuree et calculee sous des conditions P 
= 13,6 kg I cm 2 ;/= 35 Hz 


Cold-end temp Tf, K 

Net cooling capacity Qx , W 

Calculated 

Measured 

58 

0.126 

0 

69 

0.152 

0.100 

76 

0.167 

0.150 

82 

0.180 

0.200 

90 

0.196 

0.250 
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Given the refrigerator dimensions, the material phy¬ 
sical properties, the operating conditions (T L , P c u,f 
T c , T w ), and the working fluid properties, we can carry 
out the system performance analysis using the present 
model. It takes several seconds for each run on PC586. 

We carry out a performance calculation for a test 
refrigerator with design specifications listed in Table 1. 
High purity helium gas is used as the working fluid. The 
test refrigerator is designed to operate with 0.2 W cooling 
capacity at 35 Hz and 14.6 kg/cm 2 abs charge pressure. 

Table 2 shows that the measured net cooling capacity 
is very close to the calculated value from the assumption 


P ch = 14.63 kg/cm 2 abs, /= 35 Hz 



Fig. 3. Measured and calculated cooling capacity vs cold-end 
temperature. 

Fig. 3. Puissance frigorifique mesuree et calculee en fonction de 
la temperature de I’extremite froide. 


Pch~ 14.63 kg/cm 2 abs, /=35Hz 



Fig. 4. Gas bypass correction coefficient. 

Fig. 4. Coefficient de correction de la derivation du gaz. 


of r = 1 in Eq. (15) , i.e. no gas bypass in the expansion 
space. It can be seen from Fig. 3 that the predicted value 
coincides with the measured value at the design point. 
The deviation between the measured and the calculated 
values however is larger at the off-design conditions. 
This may result from the gas bypass in the expansion 
space due to the variation of displacer clearance caused 
by the temperature change. The displacer casing of the 
test ITSR is made from a plastic material that may have 
a shrinkage different from the cylinder wall material 
(stainless steel) over different Tl- 


P ch = 1.43 MPa (14.63 kg/cm 2 abs), /= 35 Hz 



Fig. 5. Comparison of measured cooling capacity with calcu¬ 
lated results. 

Fig. 5. Comparaison de la puissance frigorifique mesuree et des 
resultats calcules. 



Expansion Space Volume V e , m 3 „ 10 s 

Fig. 6. Calculated expansion space PV diagram. 

Fig. 6. Diagramme pression-volume de la chamhre d’expansion. 
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The correction coefficient r is determined by adjusting 
its value so as to fit the calculated Q L to the measured 
Q l . Fig. 4 shows the determined value of r. An empirical 
correlation is derived for P c h = 1.43MPa (14.63kg/cm 2 ) 
and/= 35 Hz: 

r=-5.60 + ^ (43) 

The predicted net cooling capacity coincides with the test 
results, as shown in Fig. 5, if the correction coefficient 
relation, Eq. (43), was used. Figs. 6 and 7 are, respectively, 
the PV diagram and the COP of the refrigerator predicted 
from the present analysis. The correction coefficient r 
may vary with the operating frequency, in addition to 



Fig. 7. Calculated COP vs cold-end temperature. 

Fig. 7. COP calcule en fonction de la temperature de I’extremite 
froide. 



Fig. 8. Prediction of net cooling capacity variation with oper¬ 
ating frequency. 

Fig. 8. Prevision des variations de puissance frigorifique nette 
selon la frequence de fonctionnement utilisee. 


the cold-end temperature. A more general empirical 
correlation of r is also derived from the analysis of test 
data: 



Fig. 9. Prediction of COP variation with operating frequency. 

Fig. 9. Prevision de la variation du COP selon la frequence de 
fonctionnement utilisee. 


Table 3 

Comparison of measured and calculated cold-end temperature 
at various operating conditions (all with P c h = 1.43MPa) 

Tableau 3 

Comparaison des temperatures d’extremite froide mesurees et 
calculees sous diverses conditions de fonctionnement (dans tous 
les cas avec Y cfl = 1,43 MPa) 


Net cooling 

capacity 

2l.1T 

Frequency, 

Hz 

Cold-end 
temperature, K 

Calculated- 
measured, K 

Calculated 

Measured 

0 

35 

57.3 

57.9 

-0.6 


40 

54.0 

56.5 

-2.5 


45 

54.0 

56.5 

-2.5 


50 

50.7 

55.0 

-4.3 

0.10 

35 

69.8 

69.2 

0.6 


40 

64.9 

64.0 

0.9 


45 

61.7 

61.4 

0.3 


50 

59.6 

59.7 

-0.1 

0.15 

35 

76.6 

75.6 

1.0 


40 

70.6 

68.6 

2.0 


45 

66.6 

65.8 

0.8 


50 

64.1 

63.8 

0.3 

0.20 

35 

84.2 

82.0 

2.2 


40 

76.8 

74.3 

2.5 


45 

71.9 

70.4 

1.5 


50 

73.4 

72.1 

1.3 

0.25 

35 

92.8 

90.3 

2.5 


40 

83.5 

79.8 

3.7 


45 

77.4 

75.2 

2.2 


50 

73.4 

72.1 

1.3 
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x (1.589 - 3.6211 x 10“ 2 /+ 5.5041 x KT 4 / 2 ) (44) 

where T l is in K; /is in Hz. 

Table 3 shows the comparison between measured and 
calculated T f of the test refrigerator. It is shown that all 
the deviations are within ± 5 K, mostly within ± 3 K. 

7. Discussion and conclusion 

The performance calculation of an ITSR may not be 
so difficult. If we look at the performance of an ITSR 
from the point of view of system dynamics, the model¬ 
ing can be greatly simplified. It is shown in the present 
study that a transfer function can be derived to repre¬ 
sent the system dynamics behavior of an ITSR and from 
which the thermal performance can be evaluated with 
sinusoidal signal assumption. The computation on 
PC586 just takes a few seconds for a design calculation. 
The performance calculation for a test refrigerator has 
shown that the analysis is accurate at the design point. 
An empirical correlation for the correction coefficient r 
in the expansion space was also derived in the present 
study for better results at off-design conditions. A PC- 
based computer package was also developed in the pre¬ 
sent study so that engineer can easily use it for design 
purposes. The modeling can be improved furthermore 
by applying a better oscillating-flow heat transfer and 
frictional correlation for the regenerator and the gas 
passage, if they are available. It is also noticeable that 
the manufacturing quality and the material selected for 
the displacer, the regenerator, and the cylinder wall are 
crucial. Any moderate defect during manufacturing will 
cause the performance to deviate largely from the 
prediction. The correction coefficient r reflects to some 
extent on this factor and the function would depend 
largely on individual overall engineering capability. 
Finally, it should be emphasized that the present study 


has demonstrated how a linear system theory or system 
dynamics approach can be applied to the design of an 
ITSR. This approach enables us to draw a performance 
map of a specific ITSR at various operating conditions and 
study the optimum design. The performance calculation 
follows the frequency-response analysis and presents the 
performance over all operating frequencies. Figs. 8 and 
9 are the prediction of Q L and COP for the present test 
ITSR operated at P ch = 14.63 kg/cm 2 abs and T’ L = 60K. 
It is seen that an optimum operating frequency can be 
determined with respect to either the maximum net 
cooling capacity Q L or the maximum COP. If we repeated 
the same calculation of Figs. 8 and 9 for different 
P ch or 7 l, we can draw a performance map for a specific 
ITSR for design selection. 
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